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O ■ Abstract. We study a dynamical curvature instability caused by a local 

^ ■ chemical modification of a phospholipid membrane. In our experiments, a basic 

solution is microinjected close to a giant unilamellar vesicle, which induces a 
2 ' local chemical modification of some lipids in the external monolayer of the 

H ^ membrane. This modification causes a local deformation of the vesicle, which 

I , then relaxes. We present a theoretical description of this instability, taking into 

'"rt ' account both the change of the equilibrium lipid density and the change of the 

^ ' membrane spontaneous curvature induced by the chemical modification. We 

Q ' show that these two types of changes of the membrane properties yield different 

(\ ' dynamics. In contrast, it is impossible to distinguish them when studying the 

equilibrium shape of a vesicle submitted to a global modification. In our model, 

the longest relaxation timescale is related to the intermonolayer friction, which 

plays an important part when there is a change of the equilibrium density in one 

^ ' monolayer. We compare our experimental results to the predictions of our model 

^+ ' by fitting the measured time evolution of the deformation height to the solution 

|^~s ^ of our dynamical equations. We obtain a good agreement between theory and 

QQ , experiments. Our fits enable us to estimate the intermonolayer friction coefficient, 

yielding values consistent with previous measurements. 
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1. Introduction 

Much more than simple barriers, membranes play an active role in many biological 
phenomena, in particular when they bend and curve into a multitude of shapes [U [2] . 
Indeed, each membrane shape is coupled to a specific function, and the diversity and 
dynamics of these shapes are vital for cell physiology [3l |4] . The remarkable material 
properties of biological membranes, inherent in their nature of complex bi-dimensional 
visco-clastic films, play an important part in membrane deformations. 

Giant unilamellar vesicles (GUVs) are model membranes made of amphiphilic 
lipid molecules that self-assemble into closed bilayers in water. Their sizes and 
curvatures are similar to those of living cells and their lipid bilayer membrane exhibits 
the basic properties of biological membranes. These features have made them very 
attractive objects to study the physics of a number of phenomena in cellular biology. 
Although lacking membrane proteins and a cytoskeleton, they have been used as a 
minimal cell model to mimic various biological processes such as membrane budding 
and endocytosis, fusion and fission, transport phenomena across the membrane, lipid 
domain formation, etc. [3 [6l [71 IH [9l [10| . Especially, giant vesicle "micromanipulation" 
gives a unique opportunity to carry out controlled experiments on an individual vesicle, 
obtaining data directly while exposing this vesicle to mechanical, biochemical, or 
chemical perturbations [TTl [T^l [131 [SI ■ In this last case, a localized chemical gradient 
is created close to the membrane by microinjection. The response of the membrane 
to such a localized gradient can be significantly different from the one caused by a 
uniform perturbation of its environment. In particular, we will show in this paper 
that a change of the spontaneous curvature and a change of the preferred area per 
lipid lead to different dynamics in the case of a local perturbation. In contrast, when 
uniform perturbations arc studied, observing the change of the equilibrium shape 
does not enable to distinguish between these two types of changes of the membrane 
properties [15]. 

In a previous work |16| . we have reported both theoretically and experimentally 
a chemically-driven membrane shape instability observed when GUVs are submitted 
to a local pH increase. In this instability, a local chemical modification affecting some 
lipids in one monolayer of the membrane changes the preferred area per lipid, thereby 
inducing a transient local bilayer curvature. Since only one monolayer is affected, 
the lateral redistribution of the lipids is strongly slowed down by the intermonolayer 
friction. When such a chemical modification is applied to a small enough surface, it 
triggers the ejection of a tubule growing exponentially toward the chemical source. 

In this paper, we describe in more detail the first phase of this local curvature 
instability, before the tubule ejection, and we present a more complete comparison of 
our experimental results with the predictions of our model. In the experiments, the 
chemical modification of the vesicle membrane has been achieved by locally delivering 
to the membrane outer leaflet a basic solution of NaOH with a micropipette. We 
have performed a time-step chemical modiflcations ( "pulse" experiments) and we have 
measured the time evolution of the height of the deformation of the vesicle with 
respect to its initial shape. Our theoretical model has been completed by taking into 
account in our dynamical equations the spontaneous curvature change in addition 
to the equilibrium density change induced by the chemical modification. We have 
also compared these two effects. The comparison of the experimental results with 
the predictions of our model has been performed by fitting the measured deformation 
height with the solution of the dynamical equations of our model. Several "pulse" 
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experiments, done on three different GUVs, have been analyzed. The agreement 
between the experiments and the model is quite good, and our fits enable us to 
estimate the intermonolayer friction coefficient, yielding values consistent with the 
literature. 

2. Materials and Methods 

2.1. Membrane composition and giant vesicle preparation 

The following lipids were used without further purification: egg yolk L-a- 
phosphatidylcholine (EYPC), Sigma, Lyon, France; Brain L-a-phosphatidylserine 
(PS), Avanti Polar Lipids, Alabaster, AL. All others chemicals were of highest purity 
grade: HEPES, Interchim, Montlugon, France; EDTA, Sigma; NaOH, Sigma. 

Giant vesicles were formed by the liposome electroformation method [17j in 
a thermostated chamber. The particular electroformation protocol used in this 
work was the following: lipid mixture solutions were prepared in chloroform/diethyl 
ether/methanol (2:7:1) with a total lipid concentration of 1 mg/ml. All the liposome 
preparations were made with a unique lipid mixture of EYPC and PS with EYPC/PS 
90:10 mol/mol. A droplet of this lipid solution (1 fil) was deposited on each of the 
two parallel platinum wires constituting the electroformation electrodes, and dried 
under vacuum for 15 min. An AC electrical field, 10 Hz, 0.26 Vpp, was applied 
to the electrodes. Buffer solution (2 ml, pH 7.4, HEPES 0.5 mM, EDTA 0.5 mM, 
temperature 25°C) was added to the working chamber (avoiding agitation). The 
voltage was gradually increased (for more than two hours) up to 1 Vpp and maintained 
during 15 more minutes, before switching the AC field off. The GUV were then ready 
for further use. In each preparation at least 10 GUVs of diameter 50-80 /im were 
available. 

2.2. Microscopy imaging and micromanipulation 

We used a Zeiss Axiovert 200M microscope, equipped with a charged-coupled device 
camera (CoolSNAP HQ; Photometries, Tucson, AZ). The experiments were computer- 
controlled using the Metamorph software (Molecular Devices, Downington, PA). The 
morphological transformations and the dynamics of the membrane were followed by 
phase contrast microscopy. 

Tapered micropipettes for the local injection of NaOH were made from GDC-1 
borosilicate capillaries (Narishige, Tokyo, Japan), pulled on a PC-10 pipette puller 
(Narishige, Tokyo, Japan). The inner diameter of the microcapillary used for 
performing the local injections onto a GUV was 0.3 /im. For these local injections, 
a microinjection system (Eppendorf femtojet) was used. The micropipettes were 
filled with a basic solution of NaOH (1 M, pH 13). The injected volumes were 
on the order of picoliters, the injection lasted a few seconds, and the injection 
pressure was 100 hPa. The positioning of the micropipettes was controlled by a high- 
graduation micromanipulator (MWO-202; Narishige, Tokyo, Japan). The injections 
were performed at different distances from the GUV surface, taking care to avoid any 
contact with the lipid membrane. The injected solution covered about 10 % of the 
GUV surface at the time of the injection. Estimations based on visualizing the flux 
from the micropipette and taking into account the dilution of NaOH in the GUV 
formation buffer (pH 7.4) yield, at the deformation onset, a pH ranging from 8 to 9 
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3. Experiments 

Giant unilamellar vesieles in the fluid phase were formed by electroformation at 25°C 
in a buffer at pH 7.4 from a mixtm-e of EYPC and PS with EYPC/PS 90:10 mol/mol. 

The chemical modification of the membrane was achieved by locally delivering to 
the membrane outer leaflet a basic solution of NaOH (1 M, pH 13). This local increase 
of the pH should affect the head groups of the phospholipids PS and PC forming the 
membrane. Indeed, the amino group of the PS head group deprotonates at high pH, 
its intrinsic pKa being about 9.8 in vesicles constituted of a PC/PS 90:10 mol/mol 
mixture [18j . Besides, the positively charged trimethylammonium group of the PC 
head group associates with hydroxide ions at high pH, the dissociation constant K 
of this equilibrium being such that pK = 14 — pK°^, where pK°^ = 11 [TS]. Both 
reactions increase the negative charge of the lipids, which entails a local change of the 
preferred area per lipid head group in the outer leaflet. 

Figure [1] shows a typical "pulse" experiment. We try to perform a time-step 
chemical modification. For this, we quickly approach the pipette without injecting 
any solution, and then we inject the solution during a time At « 3 s before quickly 
withdrawing the pipette. One can see in Fig. [1] the vesicle before any microinjection 
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Figure 1. Time-step diemical modification of the vesicle membrane. A local 
modulation of the pH at the level of a vesicle membrane induces a smooth 
deformation of the vesicle (frames 1.3 to 3.6 s). The deformation is completely 
reversible when the NaOH delivery is stopped (frames 4.3 s to the end). Scale 
bar: 50 ^m. 

(frame s), and the micropipette during its approach (frames 0.6 and 1.3 s). A two- 
step process occurs. First, a smooth deformation of the vesicle starts to develop toward 
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the pipette, i.e., opposite to the flow (frames 0.6 to 3.6 s). At tliis point we stop the 
injection and we quickly withdraw the pipette, and the membrane deformation relaxes 
after reaching a maximum (frames 4.3 s to the end). 

The following control experiments have been carried out: (i) We have checked 
that no deformation occurs if only buffer solution is injected, (ii) In order to verify 
that the observed effects were not simply due to charge screening and/or osmotic 
effects, a local injection of salt solution (NaCl instead of NaOH) has been performed. 
The typical smooth and reversible deformation has not been observed in these control 
experiments, which shows that the pH increase is crucial in our instability. 

Fig. [2] shows the measured time evolution of the height H{t) of the deformation 
of the vesicle with respect to its initial shape for the experiment presented in Fig. [1] 
One can see easily on this figure the two-step process described previously. For more 



ffi 



20 
15 
10 


♦ 
♦ 
♦ 
♦ 


1 ' 1 ' 1 ' 1 ' 

♦ 

♦ 

♦ 
♦ 


5 



- ♦ 
♦ 






500 

400 

300' 

200; 

100 



6 
t(s) 



8 10 




12 



Figure 2. Typical example of the time evolution of the vesicle deformation, 
measured in front of the pipette, in a "pulse" experiment. Dots: time evolution 
H{t) of the deformation amplitude. Solid line: Time evolution d{t) of the distance 
of the micropipette from the electrode supporting the GUV. 

clarity, the distance of the micropipette from the electrode supporting the GUV is 
also presented. 

Several "pulse" experiments were conducted on three different GUVs. The 
complete analysis of these experiments is presented in Sec. [5l 

4. Theoretical model 



4.1. Free energy of a bilayer membrane 

The free energy per unit area / of a lipid membrane is well described on a large scale 

by 



/ 



-TO + -C 



(1) 



This free-energy density depends on the membrane curvature c, which is defined as 
the sum of the two principal curvatures on the bilayer midsurface. The constitutive 
constants ctq, k and Cq denote, respectively, the membrane tension, its bending elastic 
constant, and its spontaneous curvature. The resulting membrane free energy 



F= dAf 



(2) 
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is known as the Helfrich Hamiltonian [TH] . The integral in Eq. ([5]) is a surface integral 
over the area A of the bilayer midsurface. Note that we do not include any Gaussian 
curvature term in this free energy. Indeed, the topology of the vesicle is not affected 
by the instability we wish to study, so that the integral of the Gaussian curvature, 
and thus its contribution to F, remains constant by virtue of the Gauss-Bonnet 
theorem [20] , 

To account for the bilayer structure of the membrane, its free-energy density / 
can be written as the sum of the free-energy densities of the two monolayers, which will 
be noted /+ and /~. Since the curvature c is defined on the bilayer midsurface, it is 
the same for both monolayers. Furthermore, we assume that the two monolayers 
have the same lipid composition before the onset of the instability, so that they 
have identical tensions and bending rigidities. They also have opposite spontaneous 
curvatures, noted icg, since the lipids in the two monolayers are oriented in opposite 
directions: their hydrophilic heads are oriented towards the exterior of the bilayer, 
while their hydrophobic tails are oriented towards the interior of the bilayer. To study 
our instability, it is necessary to take into account the inhomogeneities in the lipid 
mass densities p^ defined on the bilayer midsurface in each monolayer. Defining the 
scaled densities r = {p — po)/po, where po is a reference density, which is chosen 
identical for both monolayers, we may write [21j 

/ = Y+ 4C ±— c+-(r±±ec) . (3) 

In this formula, k is the stretching elastic constant of a monolayer, which is the same 
for both monolayers as they are identical, and e denotes the distance between the 
neutral surfaces of the monolayers [22] and the midsurface of the bilayer. Indeed, the 
scaled densities in each monolayer at a distance e from the bilayer midsurface are 
r^ = r^ ± ec at first order in the small variable ec, so that if /^ is written in terms 
of the curvature and r^ = (p^ — po)/po, these two variables are decoupled. Such a 
decoupling between deformations where only the curvature is modified (bending) and 
deformations in which only the density is affected (stretching) is characteristic of the 
neutral surface [22]. We choose the sign convention for the curvature in such a way 
that a spherical vesicle has c < 0. Then the monolayer denoted "-I-" in Eq. ([3]) is the 
outer monolayer. 

The expression Q of the free-energy density of a monolayer is a general second- 
order expansion around a reference state characterized by a flat shape (c = 0) and a 
uniform density p^ = po [23] . It is valid for small deformations around this reference 
state: r = 0{e) and ec = C(e)j where e is a small nondimensional parameter 
characterizing the deformation. 

If the densities on the neutral surfaces, p^, are both equal to the reference value 
Po, summing the monolayer free-energy densities /* in ([3]) gives back the standard 
free-energy density ([T]) of the bilayer: / = /^ + /^- Note that we find Cg = 
for the spontaneous curvature of the bilayer, since we have considered two identical 
monolayers. 

4-2. Modification of the free energy due to a local pH change 

In the experiment, when the pipette expels the NaOH solution close to the GUV, it 
induces a local increase of the pH, which affects the head groups of the phospholipids 
forming the membrane, as explained in Sec. |31 As some of these head groups become 
more negatively charged during this modification, the preferred area per lipid head 
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group increases locally. Pursuing the analysis of Ref. [T^, we call (j){r,t) the fraction 
of the lipids of the external monolayer that are chemically modified by the local pH 
increase in the experiment. Since the pH on the membrane should never exceed 9, this 
fraction should remain very small (see the pKa values in Sec. [3]), so we may assume 
^ = Oie). 

Since the characteristic timescale of an acido-basic reaction is determined by the 
diffusion time of the reactants [21], we can consider that (j){r, t) follows instantaneously 
the pH field outside the vesicle. The hydroxide ions diffuse in the buffer solution with 
a diffusion coefficient Dqh- ~ 5 x 10^ /im^/s [55]. Thus, in the typical 5 s of relaxation 
of the membrane deformation, they diffuse on a length of approximatively 1.5 x 10^ //m. 
This length is larger than the width of the instability and comparable with the size of 
the GUV. Besides, the pH field is probably also affected by the advection created by 
the flux coming from the micropipette, then by its retraction and last by the movement 
of the membrane. However, the central zone of the instability remains in contact with 
the zone of highest pH. At present, we have no knowledge of this time-dependent pH 
field beyond these rough estimates. Hence, to simplify, we are going to assume that 
^(r, t) does not evolve significantly during the time of the instability. More precisely, 
we make the simplifying hypothesis that (f){r) is time-independent after the end of the 
injection {t ~ 0): 

(SHi) 0(r,t) = 0(r) fort >0. 

We also assume that the inner monolayer is not affected by the pH increase outside 
the vesicle. The permeation coefficient of hydroxide ions through a lipid membrane is 
Pqh- ~ lO^'^ — 10~^ cm/s, which yields a negligible pH increase inside the vesicle on 
the timescale of our experiments [IHl HZ] . 

A priori^ the constitutive constants of monolayer "-I-" are all affected by the 
chemical modification, i.e., they depend on <^(r). However, since we focus on small 
deformations around the flat shape and the uniform density p^ — Po, we can write to 
second order in e: 

f+ = ^ + a,^+ ^-4? + a (1 + r+) 01n0 -f Jc^ + I (cq + 5o0) c 

+ ^(r+ + ec)^ (4) 

With respect to expression ([3]), the monolayer tension and spontaneous curvature of 
monolayer "-I-" have been modified according to 

^^^+ai<^+^02+a(l-l-r+)<^ln0, (5) 

Co -5- Co + cucj) . (6) 

The change of the other constitutive constants (fc, k and e) is irrelevant as far as 
the free-energy density at second order is concerned. We have assumed that the 
constitutive constants were analytical functions of 0, apart from the d{\ + r"'')01n0 
term, which corresponds to a mixing entropy term per unit area [23] . Except from 
this additional term, expression ^ is the most general second-order expansion of the 
monolayer free-energy density (in the three small variables ec, r"*" and (f) when the 
total masses of modified and non-modified lipids remain constant [33] . 

4-3. Force density in the membrane 

The force density in each monolayer of a bilayer with lipid density and composition 
inhomogeneities described by the free-energy densities (|3])-(|4]) has been derived in 
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Ref. [53]. It reads to first order in e 

p+= -fca, (r+ + ec-^0) , (7) 

Pi ^ -k d, (r" - ec) , (8) 

Pn = ctqc — kAc — ke a (?' — ^ ) ?p^0 

^ aoc - kAc - fee A (r+ - r' - ^A - — A0 , (9) 

where p^ is the force density in monolayer "±" acting in a direction i tangential to 
the membrane, while p„ = p^ +Pn is the total normal force density in the membrane. 
In these formulas, we have defined the constants k = k + 2fce^ and cq = cq + 2a\el k. 
The symbol di denotes partial derivation in the direction i, while A is the covariant 
Laplacian operator. 

Let us discuss the terms depending on in the force densities. For a given 
composition 0(r), the equilibrium density of a plane monolayer with fixed total mass 
is obtained by minimizing the free energy per unit mass f / p for c = 0. Considering 
a fixed total mass is justified here since the timescales of our instability are much 
shorter than the fiip-flop characteristic time, which is assumed not to be significantly 
modified by the local chemical modification we study. For monolayer "+" , this gives 
at first order in e 

+ ,^, _ P^qW-PO go/2 + gl0 

where we have used the fact that a^) <^ k. Thus, the chemical modification changes 
the scaled equilibrium density of the plane monolayer "+" by the amount 

5<=^c+q(</')-<(O)-y0. (11) 

We observe that this change of the equilibrium density at c = appears in the force 
densities (O and ([9]) , which is in agreement with the simpler theory of Ref. [16] . The 
equilibrium density obtained by minimization at c = 0, i.e. for a fiat membrane, will 
be referred to as the "plane-shape equilibrium density" in the following. 

Let us now focus on the term —^kcqAcJ) in (jH]), which comes from the chemical 
modification too, but which is not related to the change of the plane-shape equilibrium 
density. The membrane is at mechanical equilibrium if the force density vanishes. It 
can be seen from ([7]) and ^ that the flat shape c = is a solution to p = at a 
given inhomogcneous (f>(r) only if cq = 0. Thus, co4> (and not co<?!>) represents the 
actual change of the spontaneous curvature of monolayer "-I-" caused by the chemical 
modification. 

4.4- Change of the spontaneous curvature and of the equilibrium density 

Changing the local spontaneous curvature of a bilayer membrane may result in shape 
or budding instabilities [IHl [Ml 130] • Alternatively, affecting locally the plane-shape 
equilibrium density non-symmetrically in each monolayer of a bilayer can also yield 
shape or budding instabilities [3T1[T6]. However, until now, no study has considered 
an asymmetric modification of both the spontaneous curvature and the plane-shape 
equilibrium density in the two monolayers. As can be seen from Eq. ^, changing the 
plane-shape equilibrium density of the lipids produces a destabilizing normal force per 
membrane unit area Spn = eaiAcj), while changing the spontaneous curvature yields 
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(2) 1 - 

Spn — —^kcqAcJ). Thus, both of these changes should induce a shape or budding 
instabihty. 

In the ease where only the equilibrium density is affected, which corresponds to 
(Ti 7^ and cq = 0, the budding should vanish as soon as the lipid density has relaxed 
to its new equilibrium value, even if the modified lipids remain in place. Indeed, if 
r~ is homogeneous, and if r~^{r) reaches r|^(r) (defined in Eq. (jlOp ). the equilibrium 
condition p„ = Pi =0 is satisfied for c = 0, which means that the flat shape is an 
equilibrium shape. On the contrary, in the case where only the spontaneous curvature 
is modified, i.e., assuming cti = and cq ^ 0, the budding persists as long as the 
modified lipids remain in place. Indeed, the equilibrium condition Pn = Pi =0 can 
be satisfied only if c(r) verifies croc(r) — KAc(r) = iKCoA(/)(r), which implies that 
c(r) 7^ if A(/)(r) ^ 0, so that the plane shape is not an equilibrium shape for a 
generic inhomogeneous 0(r). Thus, although both mechanisms should lead to a shape 
or budding instability, they are not physically equivalent. We will study the differences 
between them in more detail in Sec. 14.71 

It is interesting to determine the relative importance of the two above-mentioned 
mechanisms in a budding instability. For this, let us first compare the relative variation 
of the spontaneous curvature dco/co and the relative variation of the plane-shape 
equilibrium density dp^ / p^ induced by the chemical modification in monolayer -I-. 
As explained in Sec. [3l this chemical modification affects the head groups of some 
lipids, which become more negatively charged, so that the preferred area of the head 
group of these molecules increases. This means that our chemical modification may 
be viewed roughly as an increase of the preferred area per lipid head group. It is now 
necessary to resort to a microscopic model to express the variation of the spontaneous 
curvature and of the plane-shape equilibrium density as a function of the variation of 
the preferred area per lipid head group. 

The simplest way of doing this is to take a mere geometrical model in which each 
lipid is constituted of a head group with area aho = 'ho ^^"-^ '^^ ^ chain group with 
a smaller area Cco = ^cc situated at a fixed distance s from each other (see Fig. [3]). 
Both the head and the chain are supposed to be incompressible, and to favor close- 
packing. In this crude model, the plane-shape equilibrium density is equal to m/aho, 
where m is the mass of a lipid (see Fig. [3] (a)). The spontaneous curvature is obtained 
when both the heads and the chains are close-packed, i.e. when the area per chain is 
equal to Oco and the area per head is equal to Oho (see Fig. [3] (b)): geometry yields 
Co = 2(/iio — ^co)/('cos)- If the preferred area per lipid head group Oho, or equivalently 
the corresponding characteristic length /ho is modified, it gives Sp^q/pcq = —2dlho/lho 
and 5co/co = 2(5/ho/(^coCos). Thus, given that cqs ^ 1 (since lipids spontaneously 
organize into planar membranes and not micelles), we obtain 



Scq 1 Spcq 



(12) 



Co CqS Poq 

A less crude microscopic model which we may use to express the variation of cq 
and peq caused by a change of the preferred area per head group aho was exposed in 
Ref. [35]. In this model, the head and chain of a lipid have respective preferred areas 
aho and Oco, while in the monolayer, the actual areas of these groups are Oh and Qc- 
The elastic energy per molecule is written as 

giah, flc) = ^flho ( — ^ - 1 ) + -TT'^cO ( — ^ - 1 ) > (13) 
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(a) ^ho 




(b) 



*cO 




Figure 3. Two-dimensional illustration of our simple geometrical model: the 
lipid head groups have a characteristic size ijiOi ^nd the chains have a smaller 
characteristic size IcO ■ The lipids are viewed as solid cones favoring close-packing. 
(a) Plane shape: the area per lipid is given by a^o = '?q- (b) The spontaneous 
curvature cq of the monolayer is obtained when both the heads and the chains are 
close-packed. The conical shape of the lipids has been exaggerated for clarity: in 
reality, cqs <C 1. 



where K^^ and K^ are stretching clastic constants for each group. In Ref. [35], it has 
been shown that this free energy per molecule can be used to derive the area-difference 
elasticity model of a membrane. Following the same lines, it is possible to show that 
Eq. (J13p yields the following free energy per unit area of a monolayer: 

/±- V±^c+^(r±±ec)', (14) 



4 2 2 

which corresponds to our free-energy density Eq. (jH]) in the case where the reference 
density po is taken equal to p^q (which implies ctq = |23| ) . The calculations in 
Ref. |32j provide the following expressions of cq and peq from microscopic constants: 

(^hflco + Kcaho)iaho ~ Oco) 



Co 



Pcq 



s{Kh + Kc)ahoacO 
KhacQ + A'cOho 



(15) 
(16) 



where ao is the equilibrium area per lipid on the neutral surface of the monolayer. It 
is straightforward to express the relative variations Scq/cq and Spcq/ Pcq as a function 
of the relative variation of the preferred area per head group Saw/ aw, which yields 

Scq _ Kl,al^ + gcQhO '^Peq _ Kh + Kc 1 ^Poq 

Cq "' ' 



/Vh 



KliO-coiaw — flco) Pcq JVh 

where we have used the fact that cqs = (oho ~ aco)/ao ^ 1- 



(17) 



Cos Pcq 

The values of the elastic 
constants Kh and Kc should be of the same order, and we may hint at Kc < Kh since 
the chains are able to reorganize spatially. The result (fT2| from the crude geometrical 
model is recovered when Kh — >■ oo while Kc is finite, as it should. Since s « e, the 
results ([T^ and ([T7)) coming from the two microscopic models can both be written in 
the form 

(18) 



where the order of magnitude of the numerical coefficient a is one. 



Sco 

Co 


a 
eco 


^Poq 
Peq 
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We may now compare the destabilizing pressures 6pn and 5pn caused by 
each effect: \Spn \/\5pn \ = (^K|co|)/(2e|o'i|). The relative variation of the plane- 
shape equilibrium density in monolayer "+" induced by the chemical modification 
is 5p^„/po = (5r+ — (7i4)/k (see Eq. (fTT|)). Besides, the relative variation of the 
spontaneous curvature caused by the same chemical modification is Scq/cq ~ c.q(J}/cq. 
Using Eq. (|18p and the former expressions yields |co|/|(Ti| « a/{ek). Therefore, 
we obtain |(5p„ |/|(5p„ | = aK/(2fce^). Using the well-known orders of magnitudes 
K « 10^^*^ J, e « 1 nm and k « O.lJ/m^ [52], we find out that the two destabilizing 
pressures should have the same order of magnitude. 

Thus, in the normal force density Eq. ©, the destabilizing term Spn due to the 
change of plane-shape equilibrium density should usually be comparable to the other 

(2) 

destabilizing term 6pn which comes from the change of the spontaneous curvature. 
We are thus going to take both of these effects into account in our hydrodynamic 
description of the instability. 

4-5. Hydrodynamic equations 

As we focus on small deformations with respect to the plane shape, it is convenient to 
describe the membrane in the Monge gauge, i.e., by its height z = h{x, y) with respect 
to a reference plane, x and y being Cartesian coordinates in the reference plane. Then, 
c = V^/i-|-0(e^) for small deformations such that dih = 0{e) and didjh = 0(e) where 
i,j € {x,y}. Let us denote by v (x,j/,i) the two-dimensional velocities of the lipids 
within the monolayers. Recall that the fraction of chemically modified lipids, (f>{x,y), 
is assumed to be time-independent, being fixed by the surrounding pH field modelled 
as a static field (SHi). Since we have assumed (f> = 0{e), and since the flow is induced 
by the chemical modification, we also have v = 0{e). In the Monge gauge, the force 
densities ©-([S]) become [23] 



ptix,y)= -kd,[r++e\7^h--^^j , (19) 



72^ ^i 

p- (x, y)= -kd, {r- - eV^h) , (20) 

p,{x,y) =aaV^h-k\7'^h-keV'^{r+ -r-)-!^V'^<j>. (21) 

The fraction of the chemically modified lipids may be expanded in Fourier modes: 
(j){x^y) = ^ 0qe'('^^"'"'»*'^. Since we work at first order in e, the hydrodynamic 
equations are going to be linear, which justifies the use of a Fourier expansion. Since 
the region where the hpids are modified has a well-defined width ^ q^^ (see Sec. [5]), 
we arc going to identify 4>{x, y) with its dominant Fourier mode. This approximation 
will be discussed in Sec. [HJ Furthermore, we will consider a wavevector q parallel to 
the x axis (this can be done without any loss of generality by choosing the orientation 
of the axes appropriately). We thus make a second (and last) simplifying hypothesis: 

(SH2) One-mode approximation: (j){x,y) « (pqc'^'^'^. 

Following this hypothesis, we may also write r^{x,y,t) = r^(t)e''^, v^{x,y,t) = 
v^{t)e^'^^ and h{x,y,t) = /ig(i)e*^^ at hnear order in e. 

The dynamical equations describing the joint evolution of the membrane shape 
and of the lipid densities have been developed in another context in Refs. [5T1[33]. They 
have been applied to the study of the instability caused by a local density change in 
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the case where Cq — in Ref. 
read 

2 + 

- mq^Vg - ikq I 



In terms of the above Fourier components, they 



ikq [r+ - eq^hq - y 0^ 



(c^o? 



Rq^) 



+ eq-hq) 



27]qv„ 



- 27]qv^ — b I 
+ b{v+ ~ v~ 



= 0. 



0, 
dho 



=!,»,, -4,„-g.=0. 



dt 



■iqv, 



0. 



(22) 
(23) 
(24) 

(25) 



Equations (P^ and (P5|) are generahzed Stokes equations describing the balance of 
the forces per unit area acting tangcntially in monolayer "+" and in monolayer "— ", 
respectively. The first term in each of these equations corresponds to the viscous force 
density due to the two-dimensional flow of the lipids, 772 being the two-dimensional 
viscosity of the lipids. The second term is the density of elastic forces given by Eqs. ([TO)) 
and (PO]) . The third term corresponds to the viscous stress exerted by the flow of the 
surrounding fluid, which is set in motion by the two-dimensional flow of the lipids in 
the membrane. The viscosity of the external fluid is denoted by 77. The last term is 
the stress originating from the intermonolayer friction |33| , b being the intermonolayer 
friction coefficient. Equation ([24| describes the balance of the forces per unit area 
acting normally to the membrane. Its first three terms represent the elastic force 
density given by Eq. (ETI) . Its last term is the normal viscous stress exerted by the 
flow of the surrounding fluid, the normal velocity of which matches dh/dt on the 
membrane. Finally, Equation (j25p expresses the conservation of mass at flrst order in 
e. Typical values of the dynamical parameters, used throughout, are 772 ~ 10^^ J s/m^, 
77«10-3js/m3, and6«10®-103js/m4 [MIISS]. 

Comparing these dynamical equations with the ones of Ref. [151 shows that 
—cri(t)q/k corresponds to the scalar field e(r,i) introduced in Ref. [16] to describe 
the change of the equilibrium density, which is in agreement with Eq. (|lip . Ref. [16] 
focused on the effect of the variation of the plane-shape equilibrium density, which 
corresponds to taking cg = here. Here, both the change of the equilibrium density 
and the change of the spontaneous curvature are taken into account. 



^.6. Resolution of the hydrodynamic equations 

Let us define the (symmetric) average scaled density fq{t) 
(antisymmetric) differential scaled density fq{t) = ri' ~ r~ 
Eqs. ([^ - (PH)) thanks to Eq. (pS]) and adding them together gives 



r 4- r+ and the 
Eliminating v^ in 



dra 



kq 



dt ri2q + 2ri 

while substracting them yields 



(.-. 



k ' 



kq^ 



qhq 



V 



4ry 
keq^ 



keq^ \ 



4r/ 

kq" 
~2b 



qhq 



8v ' 



/ 



V 



gig" 
26 



where we have assumed 7729^ ^ b and rjq <C b, which is well verified for n/q 
corresponding to the width of the deformation observed in the experiments. 



(26) 



.(27) 



40 /^m 
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Equation (|26p shows that fg relaxes to its equihbrium value cri(j)g/k with a very 
short timescale r = ~{f]2q + 2r])/{kq). Indeed, with typically n/q w 40 ^m, and 
with the parameters given above, we obtain r w 0.3 /is, which can be considered 
instantaneous in our experiment. 

Equation (j27p can be solved by diagonalizing the square matrix involved. In our 
experiments, we have q ^ y^aojii. Indeed, the tension of a vesicle is superior to about 
10~^ J/m^, so \/ao/fi > 10^ m^^. In this regime, the eigenvalues of the square matrix 
in Eq. ^ are 

^,«,_and72«^. (28) 



This result can be checked rapidly by noting that, in the regime where q ^ a/cto/k, 
the coefRcient keq^ /b is much smaller than all the other coefficients in the matrix, 
so that the square matrix in Eq. ([27| may be approximated by an upper triangular 
matrix. The deformation thus evolves according to 

/ig(t) = Ae-^i* + Be-^^* + J^0,, (29) 

2(70 

where the constants A and B can be determined from the initial conditions on hq and 
fq. The term KCo(f>q / (2ao) comes from the constant solution of the inhomogeneous 
equation (P?]) , which corresponds physically to the residual deformation at equilibrium. 
We find that this residual deformation vanishes if cq = 0, i.e., in the case where only 
the plane-shape equilibrium density (and not the spontaneous curvature) is modified, 
which is consistent with the previous discussions. The characteristic times 7f^ and ^2^ 
are both much longer than r, the longest one being 7-j~ , which gives a total relaxation 
time of some seconds. This slow relaxation originates from the strong intermonolayer 
friction, which is involved in the antisymmetric changes of the monolayer densities, as 
the lipids of one monolayer move relative to the other monolayer. 

4-7. Qualitative picture of the dynamics in two limiting cases 



In Sec. 14. 4[ we have shown that it is physically different to change locally the plane- 
shape equilibrium density and to change locally the spontaneous curvature. Now that 
we have studied the dynamics of the instability, we are going to analyze the difference 
between these two effects. Although both of them occur in the generic case (see 
Sec. l4.4p . we are going to study the two limiting cases when only one of these changes 
occurs. 

Let us first introduce a schematic representation of lipids, which takes into account 
the preferred area per lipid on the neutral surface, which is (by definition) independent 
of the membrane curvature, and their preferred curvature (see Fig|4]). The preferred 
shape of a lipid is represented by a cone, superimposed on the lipid. The area per 
lipid on the neutral surface (within the monolayer) is symbolized by the upper surface 
of this cone, while its preferred curvature is symbolized by the angle of this cone. 

4-7.1. Dynamics in the case where cq = We are first going to study the dynamics of 
the relaxation in the case where Cq = 0. The interest of studying this theoretical case is 
to show how modifying the plane-shape equilibrium density and not the spontaneous 
curvature can induce a transient shape instability. 

Let us consider simple initial conditions at i = 0: fq{0) = fq{0) = hq{0) ~ 0. 
This corresponds to the idealized case where the injection of the NaOH solution is 
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Figure 4. First lipid (light grey): non- modified lipid. The length £„ represents 
the preferred diameter per lipid on the neutral surface, while the angle a quantifies 
the preferred curvature. Second and third lipids (dark grey or brown): modified 
lipids. For the second lipid, the modification affects only the preferred density 
on the neutral surface, but not the preferred curvature. It is the contrary for the 
third lipid. 



local in time: then, following (SHi), we have 4>{r, t) = ior t < and (/)(r, t) = (j){r) 
for i > 0. At i = 0, the preferred area per lipid suddenly increases for the lipids of the 
external monolayer that are chemically modified, since we expect their extra negative 
charge to increase repulsion. Thus, these modified lipids are effectively compressed 
(see Fig. [lib)). 

Since the coupled dynamics of fq and hq is much slower than the dynamics of 
fq, we may consider that at t = 0^ , the equilibrium state fg(O^) = aicjyqlk has been 
reached, while fq{Q'^) = still holds, so ?'^(0+) = ^cri(f>q/k. This means that after 
an infinitesimal time, half the compression of the lipids of the external monolayer is 
relaxed. However, the flow of the lipids of the external monolayer has dragged the 
lipids of the inner monolayer, because of the intermonolayer friction, thus inducing a 
dilation equal to the compression in the external monolayer (see Fig. [5jc)). 

After this very fast first step, the time evolution of the membrane deformation is 
given by 

-72* _ p-flt 



'^it) 



0, 



eqai 



(30) 



4r] 72 - 71 

which corresponds to Eq. ([^ in the case where cq = and with the above initial 
conditions at t = O"*". The deformation appears and increases until a time imax when 
it reaches a maximum [16], before decaying exponentially with timescale 7^"^. At the 



same time, the differential density rq{t) 



decreases with timescale ^^ (see 



Fig. EKd)). The membrane curving and one monolayer sliding with respect to the 
other are two distinct responses to the discrepancy between the equilibrium densities 
of the two monolayers. Here, both of these phenomena occur in the transient regime, 
but the final state (i — >■ cxd) corresponds to a non-deformed membrane (since Cq = 0): 
the discrepancy is finally solved by the relative sliding of the monolayers, which is the 
slowest process (sec Fig. [5ljc)). 



4-7.2. Dynamics in the case where ci = Let us now discuss the opposite case, where 
only the spontaneous curvature is modified. We take the same initial conditions as in 
the previous section: fq{0) = ^^(0) = hq{0) =0. At i = 0, the preferred curvature per 
lipid suddenly increases for the lipids of the external monolayer that are chemically 
modified (see Fig.[6l^b)). Here, solving Eq. ((26|) shows that fq remains equal to zero. If 
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Figure 5. Qualitative description of the dynamics of tlic instability if cq = 0. 
The lipids that become chemically modified at t = are represented in dark grey 
(brown). The cones are dark (red) for the lipids experiencing compression or 
dilation with respect to this preferred area and in light grey (green) otherwise. 
The intensity of the grey (red) colour represents the degree of the compression 
or dilation. Here, the modification affects the preferred density on the neutral 
surface, but not the preferred curvature. 



the usual approximation (|28p is used for the eigenvalues 71 and 72 , the time evolution 
of the membrane deformation is given by 

/i,W = 0,^(1-6-^=*), (31) 

which corresponds to Eq. (P^ in the case where cti = and with the above initial 
conditions. The deformation increases exponentially with timescale 7^^ towards a 
deformed final state, where hg = KCo0g/(2cro) (see Fig. [Bl^c)). 

Note that the long timescale 'y^ involving the intermonolayer friction does not 
appear in Eq. pip. Indeed, as the plane-shape equilibrium density is not modified 
here, the lipids should not have to slide with respect to the other monolayer. In fact, 
since the membrane is curved in the final state, a small sliding is necessary as the 
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Figure 6. Qualitative description of the dynamics of the instabihty if cti = 0. 
Here, the modification affects the preferred curvature, but not the preferred 
density on the neutral surface. 



equilibrium density on the membrane midsurface is curvature-dependent, contrary to 
the equilibrium density on the neutral surface of each monolayer. This small sliding 
is visible on Fig. [DJc). This effect can be accounted for by solving our full dynamical 
equations (|27p . without using the approximation (j28p for the eigenvalues. Indeed, 
solving these equations for typical values of the parameters yields a small increase of 
the differential density fq(t) = r'^ — i"Z with timescale 7j7^, and a small contribution 
to the deformation hg with timescale 7j~ . It has been checked that this contribution 
to the deformation is negligible. 

Note that in real experimental conditions, (SHi) cannot hold indefinitely. In our 
experiment, this is mainly due to the diffusion of the hydroxide ions (see Sec. l4.2|) . Even 
if the modification of the lipids was irreversible, the modified lipids would diffuse in 
the monolayer, so the deformation would relax anyway. However, the corresponding 
relaxation timescale would be much longer than the previous ones given the small 
diffusion coefficient of the lipids in a membrane -Dup ~ 1 /im^/s [36]: the modified 
lipids would take about L^/Dup w 10'' s to diffuse on a length L k, 10^ ^m. 



Jf..l.3. Comparison of the two types of changes of the membrane properties Solving 
the dynamical equations in the two limiting cases has confirmed that a local change of 
the plane-shape equilibrium density and a local change of the spontaneous curvature 
yield different dynamics. Indeed, for the length scales we consider, when only the 
plane-shape equilibrium density is modified, the deformation increases rapidly before 
relaxing towards a non-deformed shape with a long timescale ^i (typically of several 
seconds). When only the spontaneous curvature is affected, the deformation increases 
towards a stationary curved state with a short timescale 7^ (typically shorter than 
1 s) before relaxing slowly with a timescale of about one hour. 
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In contrast, when a global modification of the environment of a vesicle is 
considered, studying the change of its equilibrium shape does not allow to distinguish 
between a change of the spontaneous curvature and a change of the preferred area per 
lipid |15j . As a matter of fact, the equilibrium shape of a vesicle with fixed volume 
is fully determined within the area- difference elasticity (ADE) model by the value of 
the combined quantity 



Aao = Aao + -eg , (32) 

a 

where Aao is the nondimensionalized preferred area difference between the two 

monolayers, a a nondimensional number involving the elastic constants of the 

membrane and Cg denotes the spontaneous curvature of the bilayer [35]. The 

shape variations observed when such global modifications are performed have been 

interpreted as coming from a change of the spontaneous curvature Cg, under the 

assumption that the preferred area per lipid was not modified [151 EH |38] . 

Thus, the dynamical study of our instability induced by a local chemical 

modification of the lipids should allow to distinguish two types of changes of the 

membrane properties, which cannot be distinguished for a static global modification. 

Since both of them should be involved generically (see Sec. l4.4p . studying the dynamics 

may enable to determine their relative importance. However, this is difficult to do 

precisely in the present experiments because of the diffusion of the hydroxide ions (see 

Sec. 



5. Comparison between the experiments and the model 

5.1. Fits of the experimental results 

In the previous section, we have exposed a theoretical model for the curvature 
instability observed in our experiments. In order to compare the experimental results 
to the predictions of this model, we are going to fit the deformation height measured 
during a "pulse" experiment with Eq. ((29)) . This equation corresponds to the general 
solution of the dynamical equations of our model, and it should describe the time 
evolution of the height H{t) of the deformation of the vesicle with respect to its initial 
shape. 

We present the analysis of several "pulse" experiments conducted on three 
different GUVs (numbered 1, 2 and 3 in the following). For each experiment, we 
take as an initial condition the last point where the pipette is present. If we call 
Hq the height of the deformation at this time, which will be referred to as i = 
as in our theoretical analysis, this initial condition can be written H{0) = Hq. The 
experimental data is thus fitted with the formula 

H{t) = {Hn-C- B)e-^<^' + Be^'^^* + C , (33) 

which corresponds to Eq. ([29| with the above initial condition. There are thus four 
free parameters in our fits: _B, C, 71 and 72. Besides, since we expect that the diffusion 
of the hydroxide ions is no longer negligible after a few seconds (see Sec. l4.2|) . we have 
carried out the fits on identical intervals of duration 5 s for all the experiments, starting 
from the end of the injection (t = 0). Thus, we eliminate the longer-term evolution 
which should no longer be well-described by our model. 

The best fits of the experimental data to Eq. ([55| are shown in Fig. [7] for one 
typical "pulse" experiment for each of the three vesicles studied. Fig. [7] shows a good 
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Figure 7. Typical examples of fits of the experimental results. Dots: 
experimental data; line: best fit to Eq. Il33|l . The dots surrounded at t = 
correspond to the initial conditions. The lines at t = 5 s represent the upper bound 
of the time interval where the fit is carried out. The insets are superposed pictures 
of the initial shape of the GUV and of its most deformed shape in each experiment. 
The definition of H{t) is recalled in these insets, (a) GUV 1 - Experiment number 
5. (b) GUV 2 - Experiment number 8. (c) GUV 3 - Experiment number 15. 



agreement between the experimental results and their best fit to Eq. ([33|) . Thus, this 
formula describes well the experimental results, which is in favor of our theoretical 
model. 

As can be seen on Fig. [71 the residual deformation at equilibrium is often close to 
zero. This residual deformation, which corresponds to the constant term KCo(/)g/(2cro) 
in Eq. (j29p , is due to the change of the spontaneous curvature induced by the chemical 
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modification (see Sec. l4.4p . Thus, if the diffusion of the hydroxide ions was much slower 
than the timescales of our instabihty, observing a negligible residual deformation would 
suggest that the change of the equilibrium density is predominant over the change of 
the spontaneous curvature. However, one cannot use this argument here, because 
the modified lipids may recover their initial state by reacquiring their protons as the 
hydroxide ions diffuse away in the solution (see Sec. l4.2p . 

The best fits of the experimental results to Eq. ([29|) provide us with an estimate of 



the time constants 7f ^ and ^2^, from which we may extract the values of the vesicle 
constitutive constants b and cto, using Eq. (|28p . However, the wavevector g, which 
comes from our single-mode description of the instability (see Sec. 14. 5p . is involved in 
the expressions of 71 and 72 (see Eq. p8))). To calculate the characteristic wavevector 
q involved in each experiment, we will measure the width W of the instability, from 
which we will deduce an estimate of q through q w n/W. 

5.2. Estimation of the width of the instability 

In order to estimate the characteristic width of the instability, we have measured the 
width at mid-height of the deformed zone of the vesicle. In practice, what we call 
"deformed zone" is the zone where the height of the vesicle is larger in the deformed 
state than in the initial state. 

We have digitized the profile of the vesicle before the approach of the pipette 
and in its most deformed state, i.e., when H is maximal, which occurs just after the 
pipette is withdrawn. Fig. Elja) shows the superposed pictures of a GUV in these 
two states during a typical "pulse" experiment. The digitized profiles at these two 
times have then been fitted to polynomials z ~ P{x) (see Fig. [Sjb)). The two fitting 
polynomials have then been substracted to get the deformation Az{x), from which it 
is straightforward to measure the width at mid-height of the deformed zone W (see 
Fig. He)). 

This estimation of W has been carried out on all of our "pulse" experiments. 
We thus have a specific estimate of the characteristic wavevector q w tt/W for each 
of these experiments. This value can be used to extract an estimate of the vesicle 
physical constants from the best fit of each experiment. 

5. 3. Measurement of the intermonolayer friction coefficient b and of the vesicle 
tension ctq 



For each "pulse" experiment, the best fit of the experimental data to Eq. ([33]), carried 
out following the method of Sec. 15.11 provides an estimate of 71 and 72 (see Table [J). 
Using the value of q obtained for this experiment as explained in Sec. 15.21 (see Table[T]), 
we may deduce the intermonolayer friction coefficient b and the vesicle tension cfq from 
7i and 72 thanks to Eq. (|28p . 

Fig. [3] shows the estimates of b obtained from 16 different "pulse" experiments 
carried out on three different GUVs. The error bars correspond to the uncertainty 
on b assuming that our model describes the experiment well. This uncertainty has 
several different origins. First, recall that we fit our experimental data on a 5s time 
interval (see Sec. 15. ip . which is a somewhat arbitrary choice. In practice, the best 
fit, and therefore the estimated values of 71 (and 72), depend on the time interval 
over which the fit is carried out. This fact can be explained by the diffusion of the 
hydroxide ions, by the change in the equilibrium state which is sometimes observed 



Dynamical membrane curvature instability controlled by intermonolayer friction 20 




40 60 
X (/im) 






15 


(c) 


/ 


^\ 


10 




/ 


W \ 


5 


/ 


\ 


n 


/ 


/ 


\ 




_^ ^--^ 



20 



40 60 

X (/um) 



80 100 



Figure 8. (a) Superposed pictures of a GUV before tlie deformation and in the 
most deformed state (GUV 2, experiment number 6). (b) Polynomial fits of the 
two vesicle shapes in (a) . (c) Plot of the difference between the deformed and the 
initial shape of the GUV. This difference is calculated by substracting the two 
fits in (b). It enables to determine the width at mid- height of the deformed zone, 
noted W. 



(see Sec. 15. 1|) . and by the fact that the width of the deformation is in fact time- 
dependent (see Sec. H]). Wc have thus varied the upper bound of this time interval 
from about 3 to 7 seconds, and taken the extremal values thus obtained for 71 (and 72) 
as the bounds of the uncertainty interval over 71 (and 72). This uncertainty interval is 
generally larger than the one estimated by the fitting software for a given time interval. 
Besides, estimating b (and ctq) also relies on our measurement of W. To determine 
the uncertainty regarding W, the measurement of W described in Sec. 15.21 has been 
carried out twice for each experiment, using the two most deformed states. All these 
factors have been taken into account in the error bars in Fig. [9l the dominating one 
coming from the choice of the time interval. 

The values we find for b, i.e., 6 = 2 — 8 x 10^ J.s/m** (see Fig. |9]), are in good 
agreement with the literature [34l El]. Besides, we can see on Fig. |9]that the values 
of b extracted from the different "pulse" experiments for each one of the three GUVs 
studied are compatible given the error bars. However, the uncertainty on b is quite 
large, as shown by the error bars in Fig. [9l Besides, our results seem to indicate that 
GUV 1 has a somewhat smaller intermonolayer coefficient friction b than GUV 2, which 
was not expected, since all of our vesicles should have the same lipid composition, and 
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Table 1 . Values of the fitting parameters 71 and 72 , and of the wavevector q 
estimated as explained in Sec. 15.21 for each "pulse" experiment. The experiment 
numbers are the same as on Fig. (9] and the thick horizontal lines indicate a change 
of vesicle. The intermonolayer friction coefficient b and the vesicle tension ctq can 
be deduced from these values (see Fig. [9] for b). 
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Figure 9. Intermonolayer friction coefficient estimated from the best fit of the 
experimental data to Eq. I I33I1 . for each "pulse" experiment. 



b should only depend on this composition. However, the existence of small composition 
variations in the experiments cannot be totally excluded. 

While the intermonolayer friction coefficient b has been deduced from 71. which 
gives the longest timescale of the relaxation (a few seconds), the vesicle tension (Jq has 
to be extracted from 72, which corresponds to a shorter timescale. The fits give values 
of 72 in the range 2 — 10s~^, so that the timescale 7^^ is about 0.1 to 0.5 s (see Table[T]). 
This is very short given that an experimental point is measured every 0.3 s. Besides, 
studying Eq. ([50)1 as a function of time shows that the term e~'^^* is of significant 
importance only before the maximum of H{t) is reached, which confirms that data 
at very small times after the end of the injection would be needed to determine 72 
well. Thus, it is impossible to extract precise values of ctq from our experiments. 
Even if more data points were available at short times, one might argue that the 
determination of ctq would still be imprecise because at short times, the pipette is still 
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being withdrawn, and the resulting flow might affect the deformation. 

We have nevertheless calculated estimates of ctq from each of our fits in the same 
way as b was determined. We have obtained tensions in the range 0.6 — 6 x 10~^ J/m^, 
without notable difference between the three vesicles (the dispersion among different 
experiments on a single vesicle is similar to the one among all the experiments). The 
order of magnitude of ao is correct, although no precise values can be extracted as 
mentioned above. It corresponds to floppy vesicles, which is in agreement with their 
qualitative observation. 

Thus, fitting our experimental results provides a satisfying order of magnitude for 
b and ao, which is an argument in favor of our model. The measurements of b are more 
precise than the ones of ctq, due to the short timescale ctq is involved in. Besides, the 
observation of a relaxation with a timescale of a few seconds which is well-described 
by a term in e"''^* shows that the change of the plane-shape equilibrium density is 
important in our instability (sec Scc l4.7p . 

6. Discussion 

Let us now discuss the possible improvements of our work. First of all, given that the 
timescale of the diffusion of the hydroxide ions is comparable to the timescale of the 
relaxation of the deformation (see Sec. l4.2p . it would be very interesting to visualize the 
time-dependent pH field during the experiment, and even better to visualize directly 
the modified lipids thanks to specific fluorescent markers. This would settle whether 
(SHi) is verified or not, and for how long. Besides, the model could be improved by 
taking into account the diffusion of the hydroxide ions, so that (SHi) would no longer 
be necessary. However, it would be necessary to make strong assumptions on the way 
the hydroxide ions arc injected and diffuse. 

Another point of our model that could be improved lies in (SH2): in practice, the 
deformation of the vesicle is not single-mode, so it would be an improvement to study 
the evolution of an initial deformation with a complete Fourier expansion. This is work 
in progress, but given the uncertainty on the spatio-temporal dependence of (f> in the 
present experiments, it is difficult to gain further insight, and at the moment, the one- 
mode approximation is the best we can do. Besides, we have used the initial value of 
the width of the deformation W to calculate the characteristic wavevector q. However, 
in reality, this width can change in time during the instability, due to dispersion, which 
does not occur for our theoretical single-mode deformation. In order to have an idea 
of the importance of this phenomenon, we have digitized the vesicle shape for all the 
experimental points in one "pulse" experiment, and calculated the width W for each 
of these points, using the method in Sec. 15.21 The corresponding results are presented 
in Fig. [TOl Wc can see that W is not constant during the instability. Its variation is 
about 90% during the five second interval used for our fits. Nevertheless, it remains 
nearly constant just after the end of the injection before increasing: using the initial 
value of W was thus the best thing we could do while keeping a constant W. 

In order to determine ctq from these experiments, it would be useful to have more 
points at short times after the end of the injection. Another further experimental 
improvement would be to use larger vesicles, and to make the injection even more 
local than in the present work, so as to be closer to the assumptions of the model. 
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Figure 10. (a) Time evolution of tiie vesicle shape during a typical "pulse" 
experiment (GUV 2, experiment number 6) after the pipette is withdrawn (bold 
dashed line: initial shape), (b) Width W of the deformation as a function of time, 
extracted from the data in (a). 



7. Conclusion 

In this paper, we have reported experiments in which a local curvature instability is 
observed when GUVs are submitted to a local pH increase. 

A theoretical model of the dynamics of the instability has been described. We 
have shown that the chemical modification of the lipids resulted in a change of the 
spontaneous curvature and in a change of the plane-shape equilibrium density. Both 
of these effects have been taken into account and compared. In our description of the 
dynamics of the instability, the intermonolayer friction plays a crucial part: it gives 
the longest timescale of the instability, as the monolayers are unable to slide with 
respect to each other on short timescales. 

Our model has been compared to the experiments by fitting the experimental 
data to the theoretical formula describing the height of the deformation as a function 
of time. The agreement between the experiments and the model is quite good. The 
intermonolayer friction coefficient can be extracted from these fits, yielding values 
consistent with the literature. Finally, we have discussed possible further developments 
of our work. 

Studying the dynamics of our instability caused by a local chemical modification 
of the lipids enables to distinguish between a change of the plane-shape equilibrium 
density and a change of the spontaneous curvature. This is impossible in the 
case of a global modification of the environment of a vesicle. The fact that we 
observe a relaxation dynamics which is well-described by a mechanism involving 
the intermonolayer friction indicates that the change of the plane-shape equilibrium 
density is important in our instability. 
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